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: ABSTRACT 

The theory of the three-dimensional rotational flow of an in- 
conpressible and inviscid fluid through an axiel turbomachine is 
described and the hydrodynanical equetions are simplified by consid- 
ering an infinite number of blsdes in esch row. The forces of the 
blades on the fluid ere treated as. non-conservetive body forces dis- 
tributed uniformly about the axis. 

Formulation of the mathemetical problem leads to one non-linear 
pertial differential equetion and tvo integral equations for the 
three velocity components. A linearized solution of these simultan- 
eous equctions for any prescribed blsede loading is based on the 
considerstion thet the vorticity sener-ted by the blades is trans- 
corted downstreem by the mean axial velocity. An iteration process 
which leads to solutions of greater accuracy ie developed by consid- 
ering for each iteration that the vorticity is transported by the 
velocities found by the previous iter-tion. 

she Bessel's functions wich occur in the Green's function solu- 
tion are replaced by thsir asymptotic values and the infinite series 
is gummed to express the solution in ciosed form. The iter:tion 
process is then adapted to mechanical cslculations by dividing the 
region of vorticity into smsli rings of rectanguler cross-section and 
detormining the influence on the velocity of a unit change of vortic- 
ity in e ch of these rinses. Once this influence is esteblished it 
is relatively easy to calculate the velocities in any axial flow 


m-chine with any :re3cribed blede loadins. 
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TH TROQUC TION 


The development of turbomachinery until recent years has resulted 
primarily from empirical methods, both in design and analysis, Recently 
theoretical methods have led to the advancement of axial flow compressors 
and turbines even though, because of the complicated mathematical problen, 
many simplifying assumptions were necessarily made and solutions were 
ap  roxinate,. 

In analogy with the three dimensionel wing theory, problems of 
flow in turbomachines may be classified as: 

(1) The Direot Problem: The direct problem of calculating the flow 
in turbomachinery is that of determining the velocity field, 
the blade forces, and the distribution of energy in the fluid 
when the blade shape, the blade speed, end the appropriate 
boundary conditions ere prescribed. 

(2) Zhe Inverse Problem: The inverse problem of calculating the 
flow in turbomschinery is that of determining the velocity 
field, the blaie shape, and the distritution of enery in the 
fluid when the blade loading, the blade speed, and the boundary 
conditions are prescribed. 

(3) The Intermediate Problem: An intermediate problem which ap- 
pears to be of interest is that of determining the velocity 
field, the magnitude of the remaining forces, end the distrib- 
ution of energy in the fluid when one blade force is prescribed 


and the blade shape is partially prescribed. 
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The direct problem arises when it is desired to investigate a 
given machine “off the design point". 

The inverse problem may be solved for the initial design but 
usually, because of the structural limitation on blade shepe (see 
appendix), it will be more prectical to formulate and solve the inter- 
mediate problem. 

In order to simplify the mathematical problem the fluid will be 
assuned to be inviscid and incompressible and the blade forces will 
be treated as body forcps, uniformly distributed through the fluid, 

60 that the flow is symmetrical about the axis of rotation and the 
vorticity is no longer shed in sheets behind each blade but is contine 
uously distributed over the region downstream of the weds. nowt. 

An axial flow machine for which the inner and outer boundaries con- 
sist of concentric cireular cylinders extending to infinity in the di- 
rection of the flow will bo considered. 

The mathematical problec: is formulated by considering the time 
rate of change, alonz a streansurfece, of the tancrential component of 
the vorticity vector. The diffic ity of this problem Lies in the 
solution of the non-linear partial differential equations that describe 
rotational fluid motion. In order to overcome this difficulty a method 
of iteration ile developed whereby solutions of any required degree of 
accuracy may be obtained. 

The first step of this iteration process provides a linearized 
solution (cf. Merble-) based on the assumption, analogous to the 
Prandtl three dimensional wing theory, that the vorticity is "trans- 


ported" downstream by the mean axial velocity and is not influenced 


~S= 
by its own induced velocities, This linearized solution is not as ac- 
curate as might be expected from the above analogy, one reason being 
that the vorticity is shed in three dimensional space instead of in 
a two dimensional sheet so that the induced velocities are likely to 
be wery large. 

For the second step of the iteration process the vorticity is 
considered trensported by the velocities found by the linearized so- 
lution of the first step. The succeeding approximations are obtained 
in the same way, in each case using the velocities of the preceding 
approximation. 

The solution is obtained by finding the appropriate Creen's func- 
tion - that is, a function G (r, 23%, ) which gives the velocity, 
consistent with the boundary conditions, induced at any point of a 
circle r, z by a unit change in tangential vorticity of a vortex ring 
at radius at an axial coordimate 6. The Bessel functions which 
arise in this solution are revlaced by their asymptotic values «nd the 
Green's function, which would occur as an infinite vem ax Bessel func- 
tions, is by this means expressed in closed forn. 

The resulting oxpression for the radial velocity, a double inte- 


gral, is modified and adopted for mechanical calculation. 
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I, Notation and Symbols 
The flow is deseribed (Fig. 1) in a cylindrical coordinate systen 
r,a> . 2, by the velocity components 
Us VW, We SGapeaverois . The sorres- 
ponding radial, tenrential, and exial 


velocity components are 





_ _ By 
a er 

a ou -_ Ow 
1.” of or 

> Jt “@.frv 
S r ae 


Coordinate System and Desig- 
nation of Yelocity and Vortisce- 
ity Components 


In addition we wlll use the followin symbols. 


a) = angular velocity of rotor 


= pressure 


absolute volocity vector 


= yorticity vector 


m1 Ol <( 3 


= force vector (force of blades on fluid) 


a 
Il. Hydrodynamic Equations 
The following equations are simplified for this case of an inviscid 


and incompressible fluid in steady, adiabatic, axially syrmetrical flow. 


Lquations of biotion: 


War von “=a — ” 
ee OG 2. = BZ 

OF r eS ee “or 7 ge 

KATY ew IV = Fe 

Power oz 


—_— = 2 
ee ee 


Continuity cquations 


m= iY += oe (3) 
io ora Ow a 
wo CUCU ee 


ie 
III, Formulation of the Mathematical Problem 
a) Development of the Basic Squations 
Because of the axial symmetry only the tangential vorticity is 

associated with the radial and axial velocities, while the radial and 
axial vorticity components are only associated with the tangential 
velocity. The tangential vorticity oconstitutes an annular vortex ring 
and by considering the deformation of this ring information concerning 
the radial and axial velocities oan be gained. It is well known (of. 
Meyer“, Marble? ) that the circulation about a deforming vortex ring is 
constant in a conservative force field. Only a radial stretching san 
oocur since the flow is symmetrical about the axis so that the constancy 
of circulation requires that in a conservative field 


he wie 
ag ‘D> 0 (4) 


In the presence of non-conservative blade forces or if radial 
vorticity is present the circwlation, and henoe the quantity a, will 
not remain constant. The law -overning its variation can be derived 
from the equations of motion. The equations of motion in vector form 


can be written 

a — 7 —_ — — 

Y Se ae VP + $W-V - = (2) 
If we taxe the curl of both sides and use the equation of continuity, 
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there results 
Y-v,R = @-7V = 7 AF 


The tengential component of this vector equation is 


wu 


2 A ic Te dv I ur IFr _ IE, 
Lap + wR += = 55-7 OR “| @ ae oe 


which, using the definitions of vorticity, can be simplified to 


Z(2)= [BE] + e- 3B] — @) 


Since the circulation is directly provortional to the quantity — , 
this aquation expresses the law governing the time change of circula- 
tion around the tangential vorticity ring. The first term on the right 
is the chenge in the axial direction of the centrifugal force, the sec~ 
ond term is the change of the redial blade force in the axial direction, 
end the third term is the chance of the axial blade force in the radial 
direction. A little thought will show that these terms, in each case, 
repregent moments tending to cause rot-tion of a particle cbout a tane- 
sential axis. It is in this manner thet a non-conservative force field 
tends to effect a change in the tangential vortex ring. 

The time derivetive on the left is taken along a streamline and 


is vritten 


(4) = egl4) 0m 3 (F] (®) 





= B= 
In order to compare the relative magnitude of the terms on the right 


we use the definition of tancential velocity and write 





ug(tl= & ~u Dw _ “ J« 4» Ku Iw 
Or\ £7 rr OFOL e are re Jz re or 
a z 
- 2 - we _ wom 
Bo\ Fr f> eee r oroe 


Remembering that for the axial flow machine the radial velocity is 
enall compared to the axial velocity and expecting from physical con- 


siderations that the velocity distribution will be smooth, it appears 


that 


“ag (F) < ~ Wee 


This inequality will be useful for the first avproximations but must be 
more closely investigated for the final apnroximations. It is conceliv- 
able that for flow that differs greatly from vortex flow the redial 
change of tancential vorticity would be of such magnitude that this 
inequality is not justified. It will however be useful to group the 


smaller terms seperately, that is, to write 
a2 / a ae ia 2 Ou e 
Fe la SE) + 22 Spt aa - ee (7) 


Using the definition of vorticity we oan write an equation for the 


radial velocity as follows 


ou la ey a/v? oF Pp J Fez 
a= 7) +f 4 = £/2(E) +2 = pe ee 


yi 
SIN 
YO 
mS 
se 
. 
“ 
ae 


(8) 
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so that the left side, if set equal to zero, would be a linear partial 
differential equation for which a solution is known. 

Since the fluid is inviscid the force of the blade on the fluid 
will act normal to the blade and hence normal to the relative velocity 


so that 


Ca. 2 eT | Fas +whfz = 0 (9) 


The radial velocity and the radial force are both anall compared to the 
other force and velocity components so that with very good approxima- 


tion 
Fz = er F 9 (10) 


Using this relation in the equations for radial velocity we have one of 


the final equations of the iteration process 
a *22 2/u Dw _ a “ey +2/ OFr_ Ou - oe (11) 
572 tile) + 5H = salA) + inl“ Fo) 1 3r ar cle a7) 


For determining a relation between the tangential velocity and the 


tancentieal farce we have the second equation of motion 


Zr = a2 (rv) rw 2(rv] = 7 Fae 
which states that the variation along a streamline of the moment of 
momentum 18 equal to the mament of the tangential farce about the axis 
of rotation. Then since we will prescribe the tangential force we 


Gan express the tangential velocity as an integral along a cylinder. 
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S$, (ro des f(a 


fal - 
= Af 2 
Oe 7 w ie (rip AA 


Tw 
faz - 


(12) 


The continuity equation, when integrated, provides an equation for 


the axial velocity in terms of the radial velocity 


3 
d 
.° Se ms Fa ae (13) 


We have developed three basis equations, Mma. ll, 12, 13, to be 
used in formulating the linearized problem, for the first approximation, 


and in the construction of an iteration procedure whereby more exact 


approximations may be obtained. 


aH 
b) The Lineerized Problem 
in order to linearize the basic equations derived in the previous 
section we shell consider axial flow, where the redial velocity is very 


snall compared to the axial velocity, 
u<<w 


end will assume that the vorticity is transmitted axially downstream 
with the mean axial velocity %, 80 that 
a 
or Ae. ) = ro = ( e) 
Furthermore,if the radial force is small, we have for caloulating the 


radial velocity the expression 


a 2 iu“ a“ 
ae BS) + SS 


+ ae) + pe ae Ko) / aa 


For calculating the linearized tangential velocity Eq. 12 becomes 


oe ee gL 


yzruL + a /4 (15) 


fa= ~-~so 


and the linearized axial velocity is, from iq. 13, 
(Z=2 
2 
seg |, +(e (16) 
fa=- oO : 
We are now able to outline a solution to the linearized problem. 
We will consider here the inverse problem and will specify the two 


force components F,, and Fg - We will prescribe F.<< Fy and F,- Fo(r.z). 


and will specify the boundary conditions anplicable to axial flow. 


=) 5% 


First the tansential velocity is determined by 


222 
f F s 
woe we, Se FY” of (15) 
9g a We 


Knowing the tancential velocity we can determine the radial velocity 








fron iq- 14 
c+ B() + 598 = Fonz) 
Z 4 [até a, ar 2 (rre)/ (1h) 
with the boundar, condition that 
a= @ OO ae eee 


ve have here an elliptic partial differential equation for the redial 
velocity. Since we know the right hand side, the unhamogeneous part, 
and have complete boundary values we can solve this equation for the 
recial velocity. This equation was solved by Marble in terms of a 
Green's function G(r,2z; «,4) which gives the velocity induced at any 
point of a circle r,z by a unit chance in tangential vorticity of a 
vortex ring of radius < at an exial coordinste &. By Marble'’s solu- 


tion the radial velocity is expressed as 





i“ = g J f(x, 8) G(r,2> x 8) Axa (16) 
where 
/ vi d 
fr«,0) = va lb srlF]t wot 2 or ky] | ; (17) 
~ = Se 
2 (ean C/, (E, °C) e a (18) 
ae |. ae 
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U, (c,r) is the linear combination of Bessel functions of order one 


Open) = /, (apy), (em - 1, CBP FY, Cenc) (19) 


and the characteristic velues En are chosen for these particular bound- 
ary conditions, Zhe quantity /, is the norm of U,(¢,r) over the in- 
terval Tr), Toe 
gg rs Ue bs i Uggerr} | (20) 
Knowing now the radial and tengential velocities we have fron 


Eq. 16 the linearized axial velocity 
fee! 
wam-f fx (ra) dp (16) 
fo Cae 


We have formulated above the linearised solution as developed by 
Marble-. It should be noted that the Groen's function as derived is 
independent of the manner of linearization provided, of course, that 
the unhomogeneous pert f(<,£) is known. If then we have another value 
of £(=x~,4), more exact then the above linearized value, we can use the 
game integral (Zo. 16) to determine a =e exact radial velocity dis- 
tribution. 

This integral, involving an infinite sum of Bessel functions, is 
extremely difficult to calculate. A solution for a particular blede 
loading was obteined by Marble’. For-the axial flow machine where the 
boundaries exclude the region noar the axis, i.e., the small values of 
the arguments of the Bessel function, we can use the asymptotic values 
of the Bessel function with good accuraqy. This makes it possible to 
aum the infinite series so that the Green's function can be expressed 


in closed form. 


“lhe 
c) Asymptotic Value of Green's Function 
An asymptotic enproximation, in closed forn,is desired for the 


Green's function 


co - €,/2£-8/ 
= a U, (Ent) VU, (én) @ (18) 
G@(r.z; a 3) . me? PF” 
ne/ 
where 
Wiénr) 2 S, (emt) ¥, (yr. )- 1 any we (19) 


a | J wh, (En 7) Iu Fen ee) - JS, (Enr.)% (ene) 


ae i Us Cente) 0” Cea (20) 
Z 


= 


and €, is found fran the boundary condition that 
ia Cone) SJ, (Ene) \ (ait) — Sent) Ge =O (21) 


ror “large* (e,r) asymptotic values of the Bessel functions are 


found in Jehnxe-2nie° D 


377 
cos (Eg k- -@ ) 
ae ee 


VF 7E,L 


Z 


‘5 


Jf CExn Te) 


Cos (en,r-Z) 


JET En (22) 


i 


Uo (Ome) 


‘ SIH 
Yc ) 2 sin (€,°- 
ee Se 
——_ /$ req? 


sin (Enl-#) 


4 
VF 7 En? 
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Yo (ene) 


16 


ea 
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the corresponding characteristics values ¢, are found from the 


SOlution of Eq. 21, with the Ressel functions replaced by their asymp- 


sotic va lucs, la(he Cee 


277 z 
in = ey : w- /, a, wen (23) 





Jstn>s these characteristic values in iq. 22 we find the canplete asymp- 


totic value of the Beasel funsticns, canpatible with the boundary con- 


ditions of the problem.. 


SS — 


rT 37 
cos 





r.-r, « 
US, (E,7/) =a Sie en, SF 
77 £77 
oS RE Ze) 
gf *. c if 
i ee oe 
TT =—- nr 
2 
aire | c. 37) (24) 
Sil? Fae, 
Y, ee) -_ 7 s T 
gue ie 
- sin (FA 
iar) = maa 
77 = ar 


The asy:aptotic value of Jy: Uy» and ),, follow casily snd the asymptotic 


value of the Green's function can be written as the infinite sun 
is | pm el 
z [Ss A oy SI OT. __ 6/777 =o eS T2-Ty 
L.-T, Lz S, 


ees 








This series can be summed ani exprossed in closed fore 


(r-r.)4(«“-F;) 


Bir 7m A= [x Wha cosh 724 Ot ee (25) 
aa 
-7, a ae 


elb- 
d) The Iteration Procedure 

Using the asymptotic value of the Green's function as derived in 
the previous section three general equations can be written for the 
iteration procedure. 

@ will use tho velocities resulting from a solution of the lin- 
eerized problem to obtain the second approximation and will obtain 
subsequent approximations using velocities of the preceding approxima- 
tion. Here we are, in effeot, assuming that the vorticity is, in each 
case, transported by the velocities of the preceding approximation. 

The radial and tangential force components are prescribed and we 
start with 

us = 0 


wv, = mean axial velocity 
If the velocities which are to be otained by the mth approximation are 
denoted by the eubeeript 7 the equations for the “th approximation vel- 


ocities are 


paz Zz 
se (Geen, - seit, 224 fp (26) 





a ns ig Ee, CO ZB 4 
f82 00 A=L, 
ie / 4,-, (4,8) G(r2;~,B)aA«af li 
2-0 os 
where : 
tppeia) = = els! +l 7 9 if Sp. na Srl 32 32] 


; cosh 7 Bae - Cos pa 
G (0% 5,8) = i/2 *& 
Z (r-7)- (0 -*%) 


cosh 7 o —-cos7 





Tee Tas 
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IV. Adaptation of the Iteration Process to 


Mechanical Calculation Methods 


The integrands of the expressions for the tangential and axial 
velocity components, mys. 26 and 28, are regular over the regions con- 
Sidered and can be integrated by any numerical method with relatively 
little difficulty. Mechanical calculations can be performed using the 
equations in thelr present for. 

The integration of the equation for the radial velocity Eq. 27 
becomes rather involved by any method, and for mechanical calcula- 
tions the form of the equation must be radically changed. An inves- 
tigation of the Green's function, G{(r,z; <./) will reveal a logar- 
ithmic singularity atol=- r,@: 2. To determine the radial velocity 
at any point r,z the integration must be carried out over the entire 
region. The logarithmic singularity of the integration represents 
the influence at r,z of an axial change of vorticity at r,z. 

It will be advantageous to perform the integration in three 
parts, one part being over the regular region, away from the singular- 
ity, and the other two parts being over a small region which includes 
the singularity. Integration over the latter region containing the 
singular point consists of evaluating a regular part and evaluating 
@ part containing a purely logarithmic singularity. With this in 


mind we write the radial velocity as the sum of three integrals. 


u = sat +o (29) 


-19- 
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eae Ce 


cosh £2 oe ee 

ff Fieye. fe [pf cebrasts 08 Fase fac 8 
cosh ‘cee CGO 

L (e704 (31) 


aa 
Pal, 2-7, 


rE Pad 


= eae F(a, 9 +f, In (2) | (e440 (32) 


Obviously the integrand for u' is finite but we need to show that 


the integrand of u'' is regular. It is easily shown that 





(mar, + FF! 

: cosh C8 - COS ————" 
lire wy 77 eA ee 
x= 7 G=2 “7 Giiie; 2-2 es. - (alien 

1a-°% 1-7 
fa =)* ~ - 
i , Ge ee) 


(33) 
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There is a singularity in this integrand, however, but it does not 


occur near « G= Pex%o The distance from r,z to the singularity, mneas= 


ured radially, is given by 








x=. = = (== Fate P Ze 
eo , 25 Fae = 
— r-f, 
sea | fe ae | He —ic | 
Teor J ae rL-T, 


From these equations this sinsularity is seen to lie outside of any 


region whose center is at r,z. The integretion for u'' therefore does 


not involve a singularity. 


The integrand for u''' obviously possesses a purely logarithmic 
singularity at A 9+ Pode 


ave 
a) Evaluation of the integrai for y!' 
Using the notation that 
2-C€ & oo Uf 246 97S zoe §=6S 
fL- SL AL TL TL 
R =e 06 6 z-% O, Za 7, pers I 


we have 


“= Tf. F(«x,8 G(L2;% A dad (35) 


The function f («,8) will be different for each iteration but the 


Green's function, G(r,2; ~,7 ) will be the same. 
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Subdivision of Flow Field 
Into Gmall Restangles 


26> 
If we divide ithe region into small reotangles as indicated in 
Fige 2 we can assuue with good accuracy that f(-<, 2 ) can be repre- 
sented in one of these sufficiently small regions n by the second 


degree surface 


Cf, = CL. * Gope = “Gi fer Cea (36) 


where four values of f(=<,/7 ) in the region n are necessary to eval- 
uate the constants ©. Suppose the four values of f( o,f ) are 
r.. * Pan Substituting these values in Dqe 36 we have four equa- 


tions from which the constants © can be evaluated. 


ae = Cy, % Cre =< 7} tr Cr3 [Par a Cra 7 he 
tne = Car + Cre Km + Car Baz + Cre Xne Anz (37) 
ns = Carr **€ Cr2 Kp; © Cop fi n3 + Cnt HKyy nj 
tra ~*~ Cap + Cm Km * Cas fBne + Cane %na fone 
from which 


toss Xp, [2a =I, for? 





Anz x<n2 Snrz2 <pn2 Gn2 
fn3 Ans [2 n3 ~n3[n3 
Ane Ane Snes =e ne 
Cm = (38) 
yo <p / Ras Xp; [Fp 
ye ~n2 re Xn2 nz 
/ Xyn3 673 Appz Cn 
/ A3ng Cae ana one 
Go = exe 


le 


isa | 4 
-= = 
wey 


fur cone@tmus <2, ~<- Cal. must be determined for each rectangle. 
a . 


it will be convenient to use f( ~,f) evaluated at the corners of the 


rectansle for this. Using the above notation we can write 


/ Ne aff (cn * One A 4+ Caz + Cna < 8B) G02; 4,8) dap (39) 
- 
A 


end can further define 


Le, , IT, G, DA a/ (2 
[[<6G, al Wf (40 
5" 2 GdxdZ 


Ina = JI), 44 Gade oF 


tnz 


Ln3 


a’ = = Feuer -"Cze ya + C, se Cie Fine/ (41) 


Here the G_, depend only on the function f(er, 4) and must be deter- 


mined for each iteration. The I, depend only on the Green's function 


4 
and can be determined once and for alle it should be noted however 
that the In, are different for each r,2; that is, we must have a cer- 
tain set of ti to use in determining the velocity at each point. 
If we denote by a the velocity at the center of themth rectangle 
and use the superscript ™” to denote the corresponding daa’ then the 
velocity u’' at the center of the mth rectangle is 

a < Ta i, + Ce ee Coatoaf (42) 


1) frr? 


oD 

The = are nunerous and are cute difficult to calculate. Howe 
ever since the Green's funttion was derived by considering only the 
fundamental equations end the boundary condition and hence is complete- 
ly independert of the preseribed blade loading or the blade shape, 
these values = can be used for any axial flow problem for which 
the hubdetip ratio is the seme as that for which the hy ere calcu- 
lated. 

Considerable simplification of the problem of calculating the 
les is posal ble because of the symmetry of the Green's function. This 
function is symmetrical about 4= 2, a fact which reduces the number of 
= by almost one-half. Furthermore it possesses another sort of sym- 
metry in 4 and 2 in that the influence of the vorticity at 6 on the 
velocity at 2 is the same as the influence of the vorticity at 2 on the 
velocity at &. This condition again reduces the number of lad by al- 
most oneehalf. This latter "reciprocal" relation can be used with the 
radial distances < and r for further simplification, but the relation 


is slightly different beosnuse of the factor [= in the Green's function. 


m 
ni 


tangles into whioh the region is divided. The rectangles should be 


Tne I_, will of course depend on the number and size of the reo- 
small in the vicinity of the region where the blade forces act and 

can be larger away from the blades. Kecordins to the linearized results 
elven by marble? it appears that in most cases it is not necessary to 
carry the integration more than about five blade span lengths (ro-r,) 
upstream and downstream of the blades, but this range may not be suf- 
ficient for the higher approximations and it is believed advisable to 


calculate the I, for a somewhat greater range. 
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The expression for uu” is best expressed in matrix form for mech- 


enical calculations. Ve define 


/ — = PY (43) 


eo 
Cer 


, 8-3, ¢ (4 4) 


co 
CTY 
3 
~~ 
SSS 
| 
Y 
a 


/ 
he le I 3¢ 


7 ise ie Es¢ 
| : ae (45) 


Using this notation we can write 
he 7 a fis rr) 
/ i = ze Gry) (4b) 
2=/2,2€ 


The four column matrices [Cag | must be evaluated for each itere- 


ations The four squaro matrices [ tat | can be determined once and for 
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0) Evaluation of the Intepral for y'' 
‘che integral u'' is the regular part of the integration over 
the small region containing the point at which the velocity is to be 


calculated 


Lee red 
. L j= / cosh WE et Pao 
a LA) © cosh ro -cos 77 (HO) LE (<-a) fdadg (32) 
wa \2 
Fe |) G9 


We have from Eq. 33 the value of the integrand et the center 


point of the region, i-.e.,at4f=r,2, as 


r-f: 
f(a, oko In Se 2 (47) 


a 
It can be shown that this point is a saddle point of the logarithmic 
function in the integrand since it represents e maximum with respect 
to ~ end a minimm with respect to /. 

It follows then that if the region is sufficiently small and if 
the function f(—,/4) is amooth we can use the value of the logar- 
ithmic term at the center as an approximate mean value (constant) so 


thet an approximate value for u'' is 


ae cre 


a oe 2), Gee a, na / f(x, 3) dx aff (48) 


4 T34, 
2-8) med 
Since this integral u'' represents the effect of only one amall 
region the above expression is sufficiently accurate unless the 
change of vorticity in this region is much greater than anywhere else 


in the flow field. 
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fe can write the integral ebove for the 7th rectangle as 


aor =f ™ oe (49) 
where 
/ /- Cos 27a 
‘TI = oe lf l2-; 0 
S Ar ar (50) 


is evaluated at the mth point and 


Z2vE rrsS 


me = [ [ f(«, 8) da dB (52) 


z-&€ 7- 
is taken over the mth reotangle. 
The J are determined once and for all as soon as the subdivisions 
of the field are decided upon. The D” must be determined for each 


iteration. 
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¢) Bvaluation of ’ the Integrel for u''* 
The intecral u''' is the singular part of the integration over 
the small region containing the point at which the velocity is being 
Calculated. 


ZrE F#S 


wal f tio 6) 2-1 In 22)? ,(ace)') da 7B (32) 
pak 


2-¢€ 


we define a new function momentarily as 


g(44r) = 46 Aas/F am 
and write the integral as 
Z¥E Ftd 
oop f G(-— 47) (nf (4 (4) Ax of (53) 


The function g(«,f,r) is regular and possesses higher order 
derivative so thet it can be expanded about the center of the rectangle 


in a Macleaurin's series in the two variables «,C. 


on =, Ayr) = ae z) + (SE), la-) ( SE) 02-4) 
ea f(FZ z) lm @ 252, (x-r) (7-2) + 


(28) (0-27) o see al 


If we substitute this expansion of g(a~,@) in the integral u'' 





and neglect terms of fourth and higher orders we have 


Le =a Zz) K, + +H, He + £38), Ks ~ 
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wnere 


Zx¥E Fred 


a = ff in| (22) Asati Ax d 
zee r4+et 
Kp =~ fedin((E3) + esl") — 


Z‘é‘G ree 


~-ff (A-2] * Inf (2 - J (ER) ] da a8 


2-e r-Jt 


The integrals involving the odd powers of (“-r) or (f -2) 
vanish because of the symmetry of the logarithmic term. The above 


integrals can be ovaluated by straishtforward integration by parts 








to clive 
se ov - aie 
= - [Sen e*4 ae —/2956€ + 4. Shan = +4 tan -/ (56) 
/ - eve" 
Ae = 2, 4S@ Ine + G24—n = Ze 5") - 
6 S* ton 5 etre an “£ (57) 
/ 2 z 
ie = - [48m Se, + See Zé -*) + 





66° tan ae ee 


‘ / (58) 
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The three constants, Kye Ape Kgs are determined once and for all 


as goon as the subdivisions of the field are decided upon. 


By straightforward differentiation we can express the deriva- 


tives of g¢(<,-,.r) in terms of f(«,) and we can write the integral 


u''' for the velocity at the rectangle m as 


‘rr? mm yy) vn 
ee = ae Ee as +E, K; 
where ye Sp a ere defincd above and where 
a .. a MM os 
2, - (9)/,,2 — AT tz) 
a 9 / i 
a i 
L. eo sa. B24 7 P= 
/ a3 
i ON /r2 
/ a*f 


oS, 


7) L241) 
Ef, "2 OPC 871 dD R®/ce 


‘a ant 


(59) 
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d) The Complete Integral for u 
ve have eValuated the integral for the radial velocity in three 
parts as given by Eqs. 42, 49, and 59. Cambining these results we 
can write @ complete expression for the radial velocity at the center 


r,z of the rectangular subdivision m. 


m ™m mm ™ | 
a a. lem [hr + Cr2 i % Cail, tne 


113 


rtm 


AG P?y™ «ws “4 (60) 


¢=42,8 

The constants C, D and = are evaluated from the unhanogeneous 
part f(~,/) of the differential equation for the radial velooity 
(Sqs- 11 and 27) and are to be determined in the manner described 
in the preceding sections. These constants must be determined for 
each iteration. 

The constants I, J and K depend only on the Green's function, 
(1.e., the boundary conditions), the hub-tip ratio, and the size of the 
subdivisions and can therefore be evaluated once end for all. They are 


used in the same manner in each iteration step. 
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¥, DISCUSSION 


A method of solving the hydrodynamic equations for the incon- 
pressible flow of an inviseid fluid through an axial flow turbomachine 
has been developed aa an iteration procedure, The next logical step 
of this approach to the problem is the evaluation of the invariant 
terms of the iteration expression for the radial velocity. Once this 
is done for the several likely hub-tip ratios it will be relatively 
easy to calculate, by mechanical means, solutions for any axial flow 
machine with any prescribed blade. loading. 

In determining the degrees of ascuracy that will be required the 
restrictions (incompressibility, etc.) which were imposed to simplify 
the mathematical problem mst be considered end evaluated. The solu- 
tion should be consistent with these restrictions and of sufficient 
accuracy to indicate the nrovcr trends of the variables. It is not 
certain that the linesrized solution meets these requirements in all 
cases. The iteration procedure can be used to determine the accuracy 
of the linearized solution and if necessary to obtain solutions of 


greater accuracy. 
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VI APPENDIX 


A Brief Comparison of the Flo d Mi: db 


The extension of the more recent axial flow solutions to apply 
to mixed flow machines would be extremely difficult and would require 
consideration of several points, not significant in axial flow, which 
are of utmost importance in mixed flow. Three essential differences 
betwoen the two problems are pointed out here in order to indicate 
which of the assumptions used in the axial flow analysis would not 
be applicable to a mixed flow analysis. 

(1) The mixed flow machine contains continuous vanes as con- 
trasted with tho rotor and stator blade rows in an axial 
flow machine. The vanes therefore cannot be "twisted" 
without introducing excessive tilt away from radiel so 
that the centrifugsl forces cause large bending moments 
in the vane and prevent operation of the machine st ex- 
tremely high speeds. <A ceonstric relation must hold bde= 
tween the relative velocity, the vane forces on the 
fluid and the shape of the vanes throughout the region 
where the vanes are present, whereas for blade rows as 
in axicl flow this relation holds only in the region of 
the narrow blede and, in fact, may be concentrated in a 
"lifting line® for a good ayoroximation, 

(2) The larger radial velocities which naturally occur in 


mixed flow prohibit the simplifying assumption used 
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for axial flow that the redial velocity is very smell 
compared to the axi-l velocity. 

(3) he boundery conditions of the mathematical problem 

of the mixed flow compressor are much more complicated 
than those for axial flow end greatly increase the dif- 
ficulty of obtaining a solution. The variables are not 
separable in this case. It should be noted too that 
the boundaries will probably be very different for each 
mixed flow problem whereas they are always essentially 
the same for axial flow. 

Several types of vanes are possible but two special vane shapes 
are likely to be of interest. 

The first of these might be called "radial vanes", These vanes 
are generated by redial lines through the axis of rotation. Here the 
redial force is zero. Radial venes are necessary in a high speed 
machine because of the high centrifugel forces. If radial velocities 
and (or) pitch angles are large ag in mixed flow tnen the angles be- 
tween the vanes -nd the hub or shroud will be acute, thus increasing 
boundary layer effects. 

The second type might be called “normal vanes". These vanes are 
generated by lines through the axis but tilted in a meridional plane 
so as to be normal to the meridional trace of the streanline at ail 
points. Hore the vane force hes two components, one acting along the 
meridional trace of the stresmsurface end one acting tangentially. 


This vane is not structurally adequate for extremely high speed rotors 
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bat does heve the sood feature thet angles between the vane end the 
hub or shroud sre right angles, a fsct wiich may ninimise boundary 
leyer effects. It is interesting to note that one force component 
eccsleretes the fluid in its path in e meridional plene and the other 
component accelerates it t-nezentially so thet it appears that no 
*“wasted" forces sre present. For the exiul flow mechine radial vanes 
and normal vsnes are about the sane. 

The selection of either radial or normal vanes will lead to 


great simplification of the mathematical problen, 





Surbomachinas 


ine tuler equetions of motion are written in vector form enc in 
cylindrical comzonent form for the isentropic flow of en inviscid 
end incompressivle fluid ected upon by non-conservative body forces. 


The flow is symuctricel about the axis. 


es — / ae 
V -o YY = - soe f- (1) 
De a ee siutti ee, is / OP 
OF 7 1! owe pag Cr 
we oY ia = ra 
yar Ol Cl oa (2) 
S _ 2. @82 ef 
“i St + W J ge a thy 


By use of 2 vector identity thece equations may be written in 


_ of the vorticity 


V. A + oF) -F (3) 


where, since the flow is isentropic end adiabatic, 


/ 


V(F) = AIP eB 
is the gradient of the total energy of the fluid. In component form 
Rq. (3) is 
A) a / aa = (F) - Fr 
wE -uf = — Fe (4) 
Ae as = 21) “Fz 
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>everal interesting relations between velocity, vorticity, total 
energy, Wane shape, etc., can be derived from these equations if we 
consider the floy between two cencantric bounding surfeces acted upon 
by venes or blades rotating with an anguler volocity ©. ‘Ye will 
e3sume on infinite number of venes (or blades) so thet the forces of 
he vane on the finid may be renresented by body forces symmetrical 
about the axis. Furthermore since the fluid is inviscid these forces 


are normal to the vane and hence to the relative velocity so that 
[V- adr} F 20 (5) 


If ve tcke the 3celer »roduct of the velocity and both sides 


of Eq. 3 these results 


oO 


7 -o(B) = VF = [0-0 F - wr Fe 
or 
Z (z] = wrk (6) 
Fron the second eqation of notion 
0 ee (7) 
#rom Sqs. 6 and 7 we gee that 
£(z) — oO 2 (rv) , where rotating vanes act 


=O) , “where no rotating vanes act 


~3 Ie 
If we consider the case where the fluid hes constant energy 
far upstream and is free cf vorticity far upstrean end make the stip- 
ulation that the total energy of the fluic is changed only iy the 
action of the moving vanes (1.e., adiabatic flow) then Zq. 8 cen be 


integrated along each streamaurface with the result that 


P P _ 
(A), 4% ie 7 upstream Pi (rv) gar upstrean 


For this case of uniform energy and no vorticity far upstrean, e 


and cory are constant far upstream so that: 


| 


7 (2 / O * upatream of the vanes 


(9) 


= p(rv) ; in regions of and downstream of the 
vanes if no forces such as stationary 
vana forees occur 


If any forces, such as those resulting from stationary Vanes or 
stators, have tencentisl components then we can only say that x = 
along the stroamsurface since in tris case wo Lirv)=o, but (rd 4 O 


Combining Aqsa. 35 end 9 we have 
x _@ ee V (ao rv} a F 
and if we mitiply both sides by the vorticity 


ee Be we f2-U0(wrr)—- _Q-F =O 


we see thet 


eo v(wrvj- Ff = (10) 
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Waine only the definition of vorticity it is easily shown that 
9 vlog Se (11) 
and hence 
o v (2) =O (12) 
and further from 2,4. 10 


LQ: a = O (13) 


Using the dofinition of vorticity in1 Zq. 9 it is seen that 


warn» @ = (Ff) (14) 
Combining “q. 14 and “%q. 3 
[V -ar}) x 2 = —-fF (15) 


For this special case we can draw the following interesting con= 


clusions: 


Upstream of the Vanes} 
1. The vorticit-: 15 zero. 
2. The totel energy of the fluid is wmiforn. 
In the Region Vhere Yene Forces Act? 
1. The vorticity vector is tangent to the vene surface. (Eqs. 5 
and 13). 
ae The vorticity vector is tangsent to a surfece of constent 


total energy. (3q. 12) 
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